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Abstract 

In this work, we introduce a degenerating PDE system with a time-depending domain for 
complete damage processes under time-varying Dirichlet boundary conditions. The evolution 
of the system is described by a doubly nonlinear differential inclusion for the damage process 
and a quasi-static balance equation for the displacement field which arc strongly nonlincarly 
coupled. 

In our proposed model, the material may completely disintegrate which is indispensable for 
a realistic modeling of damage processes in elastic materials. Complete damage theories lead 
to several mathematical problems since for instance coercivity properties of the free energy are 
lost and, therefore, several difficulties arise. 

For the introduced complete damage model, we propose a classical formulation and a corre- 
sponding suitable weak formulation in an S BV-bamework. The main aim is to prove existence 
of weak solutions for the introduced degenerating model. In addition, we show that the classical 
differential inclusion can be regained from the notion of weak solutions under certain regularity 
assumptions which is a novelty in the theory of complete damage models of this type. 

For the existence results, we had to handle the following problem: During the damage pro- 
cess it might occur that not completely damaged material regions are isolated from the Dirichlet 
boundary. In this case, the deformation field cannot be controlled in the transition from in- 
complete to complete damage. To tackle this problem, we consider the evolution process on a 
time-depending domain. In this context, two major challenges arise: Firstly, the time-dependent 
domain approach leads to jumps in the energy which have to be accounted for in the energy 
inequality of the notion of weak solutions. To handle this problem, several energy estimates are 
established by F-convergence techniques. Secondly, the time-depending domain might have bad 
smoothness properties such that Korn's inequality cannot be applied. To this end, a covering 
result for such sets with smooth compactly embedded domains has been shown. 

Key Words: complete damage, linear elasticity, elliptic-parabolic systems, energetic solution, 
weak solution, doubly nonlinear differential inclusions, existence results, rate-dependent systems. 
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1 Motivation 

From a microscopic point of vievi^, damage behavior originates from breaking atomic links in the 
material whereas a macroscopic theory may specify the damage by a scalar variable related to the 
quantity of damage. According to the latter perspective, phase-field models are quite common 
to model smooth transitions between damaged and undamaged material states. Such phase- field 
models have been mainly investigated for incomplete damage. However, for a realistic modeling of 
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damage processes in elastic materials, complete damage theories have to be considered, where the 
material can completely disintegrate. 

Mathematical works of complete models covering global-in-time existence are rare and are 
mainly focused on purely rate-independent systems [MRQ61 IBMR091 IMRZ10( IMiellj by using F- 
convergence techniques to recover energetic properties in the limit. Existence results for rate- 
dependent complete damage systems in thermoviscoelastic materials are recently published in 
[RRi2j . In contrast, much mathematical efforts have been made in understanding incomplete dam- 
age processes. Existence and uniqueness results for damage models of viscoelastic materials are 
proven in |BSS05j in the one dimensional case. Higher dimensional damage models are analytically 
investigated in |BSn4[ IMTlOl IKRZllj and, there, existence, uniqueness and regularity properties 
are shown. A coupled system describing incomplete damage, linear elasticity and phase separation 
appeared in |HKllt IHK10| . All these works are based on the gradient-of-damage model proposed 
by Fremond and Nedjar |FN96| (see also |Fre02j ) which describes damage as a result from micro- 
scopic movements in the solid. The distinction between a balance law for the microscopic forces 
and constitutive relations of the material yield a satisfactory derivation of an evolution law for the 
damage propagation from the physical point of view. In particular, the gradient of the damage 
variable enters the resulting equations and serves as a regularization term for the mathematical 
analysis. When the evolution of the damage is assumed to be uni-directional, i.e. the damage is 
irreversible, the microforce balance law becomes a differential inclusion. 

For a non-gradient approach of damage models for brittle materials we refer to |FG06t IGL091 
IBabllj . There, the damage variable z takes on two distinct values, i.e. {0, 1}, in contrast to our 
work, where z G [0, 1]. In addition, the mechanical properties are described in |FG061 IGL091 IBabllj 
differently. They choose a z-mixture of a linearly elastic strong and weak material with two different 
elasticity tensors. 

Damage modeling is an active field in the engineering community since the 1970s. For some 
recent works we refer to |Car86l IDP0941 IMie95L IMKOOL iMgTTl lFre02l [JLOSl lGUE+071 IVSLllj . 
A variational approach to fracture and crack propagation models can be found for instance in 
|BFM08[ [CFM091 IGFMlOl INeglOj ILTTT] . 

The reason why incomplete damage models are more feasible for mathematical investigations is 
that a coercivity assumption on the free energy prevents the material from a complete degeneration 
and dropping this assumption may lead to serious troubles. However, in the case of viscoelastic 
materials, the inertia terms circumvent this kind of problem in the sense that the deformation 
field still exists on the whole domain accompanied with a loss of spatial regularity (cf. |RR12j ). 
Unfortunately, this result cannot be expected in the case of quasi-static mechanical equilibrium 
(see for instance |BMR09j ). 

The main aim of this work is to introduce and analyze a complete damage model for linear elastic 
materials which are assumed to be in quasi-static equilibrium. For the analytical discussion, we 
start with an incomplete damage model which is regularized in the equation of balance of forces as 
in [ MR061 IMRZ10( IBMR091 IRR12j such that already known existence results from the incomplete 
damage regime can be applied. The basis for a weak formulation of the regularized system is a 
notion introduced in |HK11| . This notion seems well adapted for the transition to complete damage 
(see also jRR12j ). The advantage is that we can deal with low regularity solutions and we are able 
to use weakly semi-continuity arguments for the passage to the limit. In our weak formulation, 
the evolution law for the damage variable which is classically described by a differential inclusion 
becomes some kind of variational inequality combined with a total energy inequality. Nevertheless, 
we are faced with several mathematical challenges since the system highly degenerates during the 
passage. 
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The major challenge is to establish a meaningful deformation field on regions where the damage 
is not complete in the limit system. For instance, it might happen that in the limit path-connected 
components of the not completely damaged material are isolated from the Dirichlet boundary. 
In this case, the deformation variable cannot be controlled in the transition to complete damage 
with Korn's inequality. We will overcome this particular issue by formulating the problem in 
terms of a time- depending domain. The domain contains all the not completely damaged path- 
connected components of the material which still possess a part of the Dirichlet boundary. Inside 
the domain, the damage evolution is still driven by a differential inclusion. The remaining area of 
the original domain consists of completely damaged material and of material parts which are not 
completely damaged and isolated from the Dirichlet boundary. Two further complicacies arise in 
this connection. 

The first issue concerns the energy inequality. The time-depending domain approach leads to 
jumps in the energy which must be accounted for in the energy inequality of the notion of weak 
solutions as well. This issue is tackled with F-convergence techniques in order to keep track of the 
energy at jump points. 

Secondly, the time-depending domain might have very bad smoothness properties which again 
might lead to a failure of Korn's inequality. This problem is approached by proving some covering 
results for these sets with smooth domains where Korn's inequality can be applied. In this context, 
we introduce some special kind of local Sobolev spaces where we look for solutions in the limit 
system. 

This paper is structurized as follows. The next section provides an overview of the notation 
we are going to use while Section [3] developes our model first in a classical setting with enough 
smoothness properties and then in a rigorous mathematical setting by presenting a weak formulation 
with 5" -By-functions for the damage and local Sobolev functions for the deformation. It is shown 
in Theorem 3.7 that the weak notion reduces to the classical formulation when enough regularity 
is assumed. The main results, i.e. Theorem 3.10| and Theorem 3.11, are stated in Section [3. 3| while 
the proofs are carried out in the subsequent Section [4} We first solve a simplified problem in Section 



4.2, By Zorn's lemma, existence of solutions to the main problems will be proven in Section 4.3 



Complete damage models are required for a realistic description of the whole damage process. 
To the best of our knowledge, there are no global-in-time existence results of complete damage 
models of such degeneracy in the mathematical literature where the classical differential inclusion 
can be regained under some additional regularity assumptions. 
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2 Notation 



Let C M" denote a bounded Lipschitz domain, D C 5$7 a part of the boundary with li^ ^(D) > 
and T > 0. The following table provides an overview of some elementary notation used in this paper: 

C^, Ti^ n-dimensional Lebesgue and Hausdorff measure 

CI.{A;MJ^) l-times continuously differentiable function on open A C M"'+-'^ with respect 

to the spatial variable x 

dJ subdifferential of a convex function J : X — t- M U {oo}, X Banach space 

B^{A) e -neighborhood of A 

Ia characteristic function and indicator function X — )• M U {oo} with respect 

to a subset A C X 
A,mt{A),dA closure, interior and boundary of A CI 

Qt, Dt nx (0, T) and D x (0, T) 

{v = 0}, {v > 0} level and super-level set given by {x £ v{x) = 0} and {x £ v{x) > 0} 

for functions v G W^''P{Q), p > n, by employing the embedding W^'P{Q,) 

c{n) 

supp(f ) support of a function v 

Let {X, II • II) be a Banach space, / C M be an open interval and /x be a positive measure. 
The space LP(I , fi; X), 1 < p < oo, denotes the p-Bochner ^-integrable functions with values 
in X (/x-essentially bounded for p = oo, respectively). We write U'{I;X) for LP{I,C^;X). The 
subspace H'^[I; X) C L^(I; X), g S N, indicates L^-functions which are g-times weakly differentiable 
with weak derivatives in L^. Moreover, the subspace BV{I;X) C L^[I;X) consists of functions 
/ G L^{I;X) with 

essvar/(/) := inf {var/((^) \ g = f C^-a.e. in /} < +oo, 

and 

fc-i 

var/(/) := sup | ^ WfiU+i) - /(ti)|| < t2 < • ■ • < ifc with ti,t2,...tk G I for k > 2|. 

1=1 

To every / G BV{I; X), we can choose a representant (also denoted by /) with vaTi{f) < +oo. Then 
the values f{t^) := lims^f± f{s) exist for all t € / (and are independent of the representant) by 
adapting the convention /((inf /)~) := /((inf /)+) and /((sup/)"'") := /((sup/)~). The functions 
f^{t) := f{t^) and f~{t) := f{t~) are thus uniquely defined for every t £ I and do not coincide 
for at most countably many points, i.e. in the jump discontinuity set Jj. Furthermore, a regular 
measure d/ with finite variation, i.e. |d/| (/) < oo, and with values in X (called differential measure) 
can be assigned such that d/((a, b]) = f^{b) — f~^{a) for all a,b £ I with a < b, cf. |Din66| . If X 
is a finite dimensional vector space we refer to |AFP00j for a comprehensive introduction. 

If X exhibits the Radon-Nikodym property (e.g. if X is reflexive) the differenial measure decom- 
poses into d/ = f'^fj. for a (non uniquely) positive Radon measure fi and a function /^ € L^{I, /i; X) 
|MV87| . The subspace SBV{I; X) C BV{r, X) of special functions of bounded variation is defined 
as the space of functions / G BV{I; X) where the decomposition 

df = f'c' + {f+-f-)nyjf 
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for an /' G L^{I;X) exists. This function /' is called the absolutely continuous part of the 
differential measure and we also write Off. If, additionally, dff E U'{r,X), p > 1, we write 
/ G SBVP{r,X). 

For the analysis of the system given in the next chapter, it is convenient to introduce local 
Sobolev functions on shrinking sets. Let G C i}rp be a subset. The intersection of G at time 
t G [0,r], i.e. Gn (H X {t}), is denoted by G{t) := {x G {x,t) G G}. We call G shrinking if G 
is relatively open in ily and G{s) C G{t) for arbitrary < t < s < T. 

In the sequel, G C ^It denotes a shrinking set. We define the following time-dependent local 
Sobolev space: 

LtHliociG;^'') := [v.G^R'' \yte {0,T], VC/ CC G{t) open : v\u,ciQ,t) G ^^(0, t; ^^(C/; M^))}. 



Here, L^Hl^^^{G;R^) coincides with Lf^^(G;M^) (see Section 4.1.1). Lf^^(G;M^) denotes the 
classical local L^-Lebesgue space on G given by 

Ll^{G;R^) := |w : G^M^ I VF CC G open : v\v G L'^{V;R^)y 

(Note that we do not demand that G should be open.) As usual, we set L^H^ loc(^) ^t-^x ioc(^' 
At fixed time points t G (0,T), we find v{t) G H^^^{G{t);'K'^). If g > 1 we write Vf for the weak 
derivative with respect to the spatial variable as well as e{v) := ^(Vw + (S/v)'^) for its symmetric 



part. The precise definition and characterization of can be found in Proposition 4.4 

Given v G LfH^^^^{G;R^) with g > 1, we say that v = b on Dt Ci G with D (1 89. and 
T-l"'^^{D) > if for every t G (0,T) and every open set U CC G{t) with Lipschitz boundary 

v{s) = b{s) on dU Ci D in the sense of traces for a.e. s G (0,t), (2) 

is fulfilled with v := v\ux(o,t) G -^^^(0, t; H^{U; M^)). 



3 Modeling and main results 
3.1 Classical formulation 

In the sequel, let the bounded Lipschitz domain 9, C M" be the reference configuration of the 
regarded body which is also clamped at the Dirichlet boundary D C dQ with 'H"'^^{D) > 0. Even 
though more general forms are conceivable, we confine ourselves in this work to the following free 
energy density f and dissipation potential density g based on the gradient-of-damage theory: 

<^{e,z,Vz) ■.= ^\Vz\P + W{e,z), f5(i) := -ai + + /(_^^o] (^), (3) 

where e denotes the linearized strain tensor, z the damage phase-field variable, a > and (3 > 0. 
The function W represents the elastic energy density. More general free energies incorporating 
appropriate convex potentials for the damage function can be employed with the obvious modifica- 
tions but we confine the analysis to the case ([3| in order not to overburden the presentation. The 
gradient exponent p satisfies p > n. The damage variable z specifies the degree of damage at each 
reference position x G 9. in the material, i.e. z{x) = 1 stands for an undamaged and z(x) = for 
a completely damage material point whereas intermediate values represent partial damage. Fur- 
thermore, the irreversibility of the damage process (the solid can not heal itself) is ensured by the 
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indicator function in q. We assume that no volume forces are acting. The flow rule yields the 
pointwise inclusion 

OG -(l\\{Lp^yz{e,z,Vz)) + Lp^z{e,z,Vz) + dIyQ^oo){z) + diQ{dtz). (4) 

Note that the subdifferential of the indicator function -/^[o.oo) appears on the right hand side of the 
inclusion to account for the constraint z > 0. 

This paper will cover elastic energy densities of the form 

W{e,z)='^-g{z)Ce:e (5) 

with a symmetric and positive definite stiffness tensor C G /^(M"y^) and a function g G C"^([0, 1]; M"*") 
with the properties 

r]<g{z), 5(0) = (6) 

for all z £ [0, 1] and some constant rj > 0. 

We use the small strain assumption, i.e. the strain calculates as 

e = e(n):=^(Vn + (Vn)T), (7) 

where the right hand side denotes the symmetric gradient of the deformation field u. 

By neglecting inertia effects, the momentum balance equation is a quasi-static mechanical equi- 
librium and reads as 

= div{W,e{e,z)). (8) 

Note that complete damage is possible if and only if g(0) = 0. The case ^(0) > would 
describe incomplete damage which is already covered in the mathematical literature (see Section 
[T]). As mentioned in the introduction, a regularization scheme is adapted, where a regularized 
elastic energy density W^, e > 0, is used instead of W in the first instance. More precisely, is 
given by 

W'{e,z) = ^ig{z) + e)Ce:e. (9) 

In the complete damage regime e = 0, the deformation variable becomes meaningless on material 
fragment with maximal damage because the free energy density vanishes regardless of the values of 
u. Therefore, the balance laws Q and ([s]) make obviously only sense pointwise in {z > 0}. Beyond 
that, as already mentioned in the introduction, a phenomenon (in the following called material 
exclusion) might cause severe troubles. 

Suppose that at a specific time point t, a path-connected component P (relatively open in il.) 
from {x G O I z{x, t) > 0} is isolated from the Dirichlet boundary, i.e. 'H"'~^{Pr\D) = 0. In this case, 
Korn's inequality fails on P and, consequently, the deformation field Us for the regularized system 
cannot be controlled on P in the transition e — )• 0"*". To overcome this problem, path-connected 
components P of the not completely damaged area {z{t) > 0} isolated from the Dirichlet boundary, 
i.e. 'H"~^(Pn D) = 0, will be excluded from our considerations in our proposed model. On the 
one hand, this make our model accessible for a rigorous analysis and, on the other hand, for some 
applications, the detached parts might be of little interest anyway. This approach is illustrated in 
Figure 1 and motivates the definition of maximal admissible subsets. 
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Complete damage model (illustration) 



Example 1: Material exclusion 




t=ti : exclusion of an undamaged part 



t=t, : after exclusion 



Example 2: Decomposition of aF(t) 




Caption 

slnrinking set F(t) 
// Dirichiet boundary D 



boundary aF(t) 
Neumann boundary an\D 



excluded parts (with no Diriciilet boundary) 
completely damaged parts 
not completely damaged Dirichiet boundary r-^(t) 
not completely damaged Neumann boundary r2(t) 
completely damaged boundary r^{t) 



Figure 1: The first example in the illustration above shows the exclusion of an undamaged material 
part during the evolutionary process in 2D. The dark blue curve encircles the maximal admissible 
subset F{t) = 2l/5({2:(t) > 0}) of the not completely damaged area {z(t) > 0}. The second example 
below pictures the different parts of the boundary of F{t). There, the Dirichiet boundary D consists 
of two components. 

Definition 3.1 (Admissible subsets of $7) 

(i) Let F <^ 0, be a relatively open subset and 

Pf(x) := {y € F\x and y are connected by a path m F } 

for x G F. We say that F is admissible with respect to the Dirichiet boundary D if for every 
X £ F the condition 

n''-\PF{x) n D) > 

is fulfilled. Furthermore, ^d{F) denotes the maximal admissible subset of F with respect to 
D, i.e. 

2l£)(F) := 1^ {G C F I G is admissible with respect to D} . 
(a) For a relatively open subset F C Qj^, the set 21/) (F) is given by (2lj:)(F))(t) := 2l£)(F(t)). 
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In a nutshell, the evolutionary problem Q and ([s]) is considered on a time-depending domain 
(a shrinking set) which is, for any time, admissible with respect to D. The whole evolution problem 
with its initial-boundary conditions can be summarized within a classical notion in the following 
way. 

Definition 3.2 (Classical solution) A pair of functions {u, z) defined on an admissible shrinking 
set F C with 

F{t)=^n{f]Fis)) (10) 

is called a classical solution to the initial-boundary data [z^, b) if the following properties are satis- 
fied: 

(i) Regularity: 

u G el{F;W), z G C2(F;M) 

(ii) Evolution laws in F: 

G -div(|Vz|P-2vz) + W,Mu), z)-a + /Sdtz + dI^_^^o]{dtz) 
0< z 

(Hi) Initial-boundary conditions: 

z{t = 0) = z°, on F{0), 

u{t)=b{t) onTi{t) 

W^e{e{u{t)),z{t)) -v = Q onT2{t) 

z{t) =0 on T^{t) 



F{t) n D, 

F{t)n{dn\D), 

dF{t)\Fit) 



Vz(t)-i/ = onTi{t)ur2{t) 

Remark 3.3 The time- depending boundary dF{t) disjointly decomposes into Ti{t)[JT2{t)L)r3{t), 
where Ti{t) indicates the not completely damaged Dirichlet boundary, T2{t) the not completely 
damaged Neumann boundary and T3{t) the completely damaged boundary (see Figure 1). We have 
the following types of boundary conditions: 

Ti{t) — Dirichlet boundary condition for u 
Neumann boundary condition for z 

^2{t) — Neumann boundary condition for u 
Neumann boundary condition for z 

T'i^t) — degenerated boundary condition 

On the degenerated boundary, z vanishes (homogeneous Dirichlet boundary condition for z ) and, 
therefore, if we assume that e{u) can be continuously extended to Fs the stress W^e{^{u), z) vanishes 
too. 

The goal of the next section is to state a weak formulation such that existence can be proven. 
Due to the high degree of degeneracy and the non-smoothness of F, u can only be expected to be 
in some local Sobolev space on the shrinking set F introduced in ([T]). 
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3.2 Weak formulation and justification 



For a weak formulation of the system presented in Section 3.1, we will take advantage of the 
free energy £ whose density has already be given in ([s]). In contrast to |KRZ11| for incomplete 
damage models and related works, we will not use a purely energetic approach but rather a mixed 
variational/energetic formulation as presented in |HK11| . 

Definition 3.4 (Free energy) Let e G I? {VL;W^:^^) and z £ W^'P{n) be given. The associated 
free energy of the system in Definition \3.^ is given by 



£(e,z):= f -\Vz\P + W(e,z)dx, 
JqP 

whereas its e-regularization with e > (for later use) is defined as (see ([9] 

£e{e,z) : 



[ ^\Vz\P + W%e,z)dx. 
JnP 



If e is only defined on a measurable subset H d VL, i.e. e G -L^(i7; Mgy^), we use the convention 
£{e, z) := £(e, z), where e := e in H and e := in i}\H . 



We are now able to give a weak formulation of the system in an S'S^-setting (with respect to the 
damage variable). In accordance to Definition 3.2 z is extended on whole and when viewed as 
an SBV^{0, T; L^(r2))-function has a jump at t if and only if a material exclusion occurs at t. 

Definition 3.5 (Weak solution) A pair (u, z) is called a weak solution of the system given in 
Definition 3.2 with the initial-boundary data {z^ .,h) if 

(i) Regularity: 

z G L-(0,T; ^^^'^'(f])) nSW2(0,r;L2(l^)), u G L^,HI,,,{F;R^) 

with e(n) =: e G L^{F;R'^^^) where F := ^dHz' > 0}) C $7 2" is a shrinking set. 
(a) Quasi-static mechanical equilibrium: 

= [ WAe{t), z{t)) -.eiOdx (11) 

JF{t) 

for a.e. t G (0, T) and for all C G H\){Q.; M"). Furthermore, u = b on DtDF. 
(Hi) Damage variational inequality: 

\Vz{t)\P-^Vz{t)-VC + W^,{e{t),z{t))Cdx> [ {a - f3d^z{t))C dx (12) 
F{t) ' Jn 

< z{t) in n, 
> dfz{t) in n 

for a.e. t £ (0, T) and for all C, G W^'P{Cl) with C < 0. The initial value is given by z^{0) = z^ 
with < z° < 1 inn. 
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(iv) Damage jump condition: 



z~^{t) = z (t)lp'(i) in O 



(13) 



for all t G [0,r]. 
(v) Energy inequality: 



f(e(t),z(t)) + 



/*/ a\dfz\+P\dfz\^d{x,s)+ y J, 

Jo Jf(s) ^ rr;„ „ 



JF{s) 

W,e{e,z) : eidtb)d{x,s) 



seJzn(o,t] 



10 JF{s) 

for a.e. t £ (0, T), where the jump part Js satisfies < and is given by 

Js := lim essinf£:(e(i?),z(t?)) - e+ 

T-s>s- i?e{T,s) 



(14) 



(15) 



and the values £ M+ satisfy the upper energy estimate 

tt <£{e{bis) + C),z+{s)) (16) 
for all C e Hl^{{z+{s) > 0};M") with ( = on D D {z+{s) > 0}. 

C.l ensures that for all times t we have z~{t)tp(^f) S W^'^{Q) (see 



Remark 3.6 (i) Lemma 
jump condition (13)^. 

(a) Jump condition ([l3]) and the definition of F imply {z~^{t) > 0} = F{t) for all t G [0,T]. By 



the convention introduced in Definition 3.4 
£{e{t),z{t)) 



F(t)P 



Vz{t)y' + W{e{t),z{t))dx, 



which equals j ^^^^-^^^^ z{t)\P + W {e{t) , z{t)) dx for a.e. t£{0,T). 

(Hi) The jump term Js equals the energy of the excluded material parts at time point s, i.e. Js = 
£{s~)—£{s~^) (for smooth solutions on F), where 1 1— )• iS(e(t), z{t)) denotes the energy function 
along the trajectory. However, for less regular weak solutions as in Definition 3.5. the one- 
sided limits £{s~) and£{s~^) possibly do not exist. But, in any case, lim^^^- essinf^g^^ £{'&) 
clearly exists and coincides with £{s^) for smooth solutions. The value £{s^), on the other 
hand, can be avoided in a rather indirect way by using upper energy estimates. More precisely, 
it turns out that £{s~^) can be substituted by values (denoted by cf) merely satisfying (16). 
Together with equations (11)-(14), is forced to coincide with £{s^) for smooth solutions. 
This is particularly shown in the proof of the following theorem. 



Theorem 3.7 Let {u,z) be a weak solution according to Definition 3.5 We assume the regularity 
properties u S C^(Ot;M"') with u = b on Dt and z\-p G C^(-F;M). Then, {u\F,z\-p) is a classical 



solution according to Definition 3.2 
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Proof. The first equation and the last inequality in (ii) from Definition 3.2 as well as property (iii) 
follow immediately by classical integral calculus from the weak notion. 



By the monotonicity of z with respect to t (coming from > dfz and the jump condition (13)) 



and by Remark 3.6 (ii), 

F{t) = ^D{{z-{t) > 0}) = 21d ( f]{z+{s) > 0}) = 21d ( fl F{s] 



s<t 



s<t 



Therefore, condition (10) from Definition 3.2 is shown. 



Finally, we need to prove the differential inclusion in (ii) from Definition 3.2 The jump condition 

nption yi 

dfz{x,t) 



(13) and the regularity assumption yields for a.e. {x,t) G fi-r 

' dtz{x,t) if(x,t)GF, 







if (x, t) G J7t \ F, 



where dtz{x,t) is the classical time-derivative of z at {x,t). In the following, we will make use of 
this property. First, observe that by the regulartiy assumptions q := {e,z) £ SBV{0,T; X) with 



Applying the chain rule (see Corollary B.2) for the continuously Frechet-differentiable energy 



functional £ and the X-valued SBV-function q shows that iS o g is an SBV-function and 
£{q{t+))-£{q{0+))=d{£oq){{0,t]) 

{de£iq{s)),dte{s)) + {d,£iqis)),dfz{s))ds 



The two terms in the integral on the right hand side can be treated as follows. 



Taking into account z = in \ F and testing (11) with C = dtu{s) — dtb{s), yields 



{de£iq{s)),dte{s))= / W,eieHs)),zis)) : eidtu{s))dx 
Jn 



Fis) 



W4€{u{s)),zis)) : eidtu{s))dx 



W4€{u{s)),zis)) : €{dtb{s))dx. 

F(s) 

• Using the property d^z = in ^It \ 

{d,£{q{s)),d!z{s)) 

\Vz{s)\P-^Vz{s) ■ Vdfz{s) + W,Mu{s)),z{s))dfz{s) dx 



F{s) 



\Vz{s)\P-^Vz{s) ■ Vdtz{s) + W^Muis)), z{s))dtz{s) dx. 
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Putting the pieces together, we end up with 



sGJzn(o,t] 



^('Z(0+)) + f f 
Jo Jf 



JF{s) 



W4eiu),z) : e{dtb)d{x,s) 



+ [ [ \Vz\P-^Vz -Vdtz + W^^{e{u),z)dtzd{x,s) 

Jo Jf{s) 



(17) 



Note that we have <5(g(0+)) = . Indeed, passing t 0+ in ([14]) yields <S(g(0+)) < Cq . The 
">"-inequahty follows from (16) tested with C, = u{0) — 6(0). 



Therefore, ( 14 ) particularly implies 

ft 

lF{s) 



S{Q{t^))+ f ! a\dtz\+p\dtz\''d{x,s)+ V Js 

<£{q{0+))+ [ [ W4e,z):e{dtb)d(x,s) 
Jo Jf(s) 



(18) 



Integrating (12) on [0,t] with respect to time, testing it with ( = dfz < 0, applying it to (17) and 



comparing the result with the energy inequality ( 18 ) shows 

ft 



£{qit+))+ V {£{q{s-))-£{q{s+)))+ [ [ -adtz + P\dtz\^d{x,s) 

> £{q{0+)) + [ [ W,e{e{u), z) : e{dtb) d{x, s) 
Jo Jf{s) 

>£{q{t+))+ y Js+ [ [ -adtz + ^\dtzfd{x,s). (19) 
^ TT^rr^u-^ Jo Jf(s) 



sGJzn(o,t] 



Taking also (15) into account and using £{q{s )) = lim^_>g- ess inf^g(T-^) <S((7(7?)), estimate (19) 
yields 



seJz seJz 



(20) 



On the other hand, by (16), we find < £{q{s~^)) for all s £ Jz- Combining this with (20) shows 
£{q{s+)) = e+ for ah s G J^. 

Therefore, = £{q{s~)) — £{q{s^)) and (19) becomes an equality. Taking also (17) into 
account gives 

= [ [ \Vz\P-'^Vz-Vdtz + W,z{e{u),z)dtz-adtz + p\dtz\^d{x,s). 
Jo Jf{s) 



Together with the variational inequality (12) and the regularity assumptions, we obtain 



< / ( - div{\Vz{s)r^Vz{s)) + W,Mu{s)), z{s)) -a + pdtz{s) ] (C - dtz{s)) dx 

If{s) 
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for all s G (0,r) and for all C G L^{F{s)) with C < 0. This leads to 

< ( - div(|Vz|P-2vz) + W^Mu),z) - a + pOtz^ {( - dtz) 

for a.e. (x, t) £ F. By the regularity assumptions, this inequality holds pointwise in F. Therefore, 
the differential inclusion in Definition 3.2 (ii) is shown. □ 

One main goal of this work is to prove ex- 
istence of weak solutions according to Defi- 
nition 3.5 Due to the application of Zorn's 
lemma used in the global existence proof, 
analytical problems arises when infinitely 
many exclusions of material parts occur in 
arbitrary short time intervals in the "fu- 
ture", i.e. cluster points from the right of 
the jump set Jz* (denoted by C^^. in the 
following) where z* G SBV{Q,T]L'^{n)) is 
given by z*{t) := z{t)t<^^i^{z-{t)>o})- In this 
case, we are only able to prove that the 
shrinking set F is approximately given by Figure 2: An example of a shrinking set where in- 
'^d{{z~ > 0}) whereas the strain e can still finitely many exclusions during an "infinitesimal" 
be represented as the symmetric gradient of time-interval have occured. 
u in 2lz)(F). 

To be precise, we introduce the following notion. 

Definition 3.8 (Approximate weak solution) A triple (e, n, z) and a shrinking set F C Qj- is 
called an approximate weak solution with fineness t] > of the system according to Definition 3.2 
to the initial-boundary data {z^ ,b) if 

(i) Regularity: 

z e L~(0,r;Tyi'f(n))n5W2(0,r;L2(J7)), u G LlHl,,,{%D{F)-mP), 
eGL\F;K;^) 

with e = e(n) in ^£){F). 

(ii) Shrinking set properties: 

F{t) D %D{{z'{t) > 0}) for all t € [0,r], 

F{t)=^D{{z~it)>0}) for allt£[0,T]\ [j [t,t + r]), 

tec,* 

C''{Fit)\^D{{z-{t)>0})) <r] for alltG \J [t,t + ry). 

teC'* 



(Hi) Evolutionary equations: 



Properties (ii)-(v) of Definition 3.5 are satisfied. 



Remark 3.9 If an approximate weak solution {e,u,z) on F according to Definition 3.8 satisfies 
Cz* = then {u,z) is a weak solution according to Definition 3.5. 
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3.3 Main results 



Theorem 3.10 (Global-in-time existence of approximate weak solutions) 

Let b G W^''^{0,T; M")) and E W^^P{n) with < < I in Q and {2° > 0}, admissible 

with respect to D, be initial-boundary data. Furthermore, letr]>0 and W be given by ^ satisfying 
Q. Then there exists an approximate weak solution {e,u,z) with fineness rj > according to 
Definition 3 



Theorem 3.11 (Local-in-time existence of weak solutions) 

Let b G H^i'i(0,r;Pri'°°(17;M")) and z^ G W^'P{n) with < k < z° < I in n be initial-boundary 
data. Furthermore, let W be given by ([s]) satisfying Q. Then there exist a maximal value T > 
with T <T and functions u an d z defined on the time interval [0,T] such that {u,z) is a weak 



solution according to Definition 3.5 Therefore, if T < T, {u,z) cannot be extended to a weak 
solution on [0, T + e] . 



4 Proof of the main result 

4.1 Preliminaries 

4.1.1 Covering properties 

The aim in this subsection is to prove covering results for shrinking sets. 

Definition 4.1 (Fine representation) Let H 0, be a relatively open subset. We call a count- 
able family {U^} of open sets CC H a fine representation for H if for every x €z H there exist 
an open set U C M" with x £ U and an k & N such that [/ n $7 C J/^. 

Remark 4.2 Note that H dO, is not covered by {Uk}. See Figure 3 for an example. 

Lemma 4.3 Let G C Vtrp be a relatively open subset and the sequence {tm} containing T be dense 
in [0, T]. Furthermore, let {U^^ken be a fine representation for G{tm) for every m € N. Then, 
for every compact set K <^ G there exist a finite set / C N and values G N, k £ L , such that 

Proof. To every element p = {x, t) G K, we will construct a neighborhood 0p C of p in the 
subspace topology of such that there exists k,m £ N with Op n fix ^ UJJ^ x {0,tm)- Then the 
claim follows by the Heine-Borel theorem. 

Indeed, to every p = {x,t) £ K there exists a e > such that Bs{p) n Qt ^ G since G C Qrp 
is relatively open. Therefore, if i < T, (x, tm) £ G for all m G N such that t < tm < t + e. This 
implies {x,tm) £ G D {0, x {tm}) = G{tm) x {tm}- Then, we find p £ G{tm) x J with J = [0, tm)- 
In the case t = T, it holds p £ G{T) x J with J = [0, T]. Since {UJI^}k£n is a fine representation of 
G{tm), let 6 > such that Bs{x) n Q C [7™ for some A; G N. Finally, Op := {Bs{x) n H) x J is the 
required neighborhood of p. □ 
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Figure 3: The left illustration shows a fine representation for the relatively open subset H = 
(0, 1) X (0, 1] of Cl = [0, 1] X [0, 1] whereas the right picture does not show a fine representation for 



H. 



A simple consequence of Lemma 4.3 is LfH^^^^{G; 



t>N\ 



Ll^{G;R^), provided that G O 



is shrinking. Moreover, we can characterize the function space LfH^^^^{G;M ) as fohows 



Proposition 4.4 Let G C be a shrinking subset and let {tm} and {UJ^} be as in Lemma 4-3 
Furthermore, let v : G ^ be a function. 

(a) The following statements are equivalent: 

(V V € LjHl,^^{G;R^) 

(n) ^^lf/rx(0,t„0 e L^{Q,t^-H\Ul^;R'')) for all A;,m G N 



(Hi) V G Lj^Q^(G;M^) and there exists a function g G Li^^^{G; . 



t)Nxn\ 



such that 



V ■ div(C) d(x, t) 



IG 

t)Nxn\ 



g : Cd{x,t) 



(21) 



G 



for allC eC^{int{Gy,: 
If one of these conditions is satisfied we write Vv := g and e{v) := ^(Vv + Vv"^). 
(b) Assume that each C/™ has a Lipschitz boundary. Then the following statements are equivalent: 
(i) V = b on the boundary Dt H G 

(a) for every /c,m G N, condition ^ is satisfied for U = [7™ and t = tm 

Proof. 

(a) (i)^=^(ii) and (iii)=^(i) are trivial. 

(ii)^=^(iii): Let the function g : G ^ M^^" be £""'"^-a.e. defined as follows. For each 



k,m en, we set g\ujj- := g^ where g^ G L^{UJP x (0,t„);: 
^lc/™x(o,t,„)- The function g is well-defined on G H Qt since 

GnnT= U urxio,tm) 

fcjTTieN 



pNxn 



) is the weak derivative of 
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and = on C/^^ x (0,t^J n [/^^ x (^t^J for all A:i,fc2,mi,m2 G N in an £"+i-a.e. 
sense. Let t G (0,T] and [/ CC G{t) be open. By Lemma 4.3, U x (0, t) can be covered 

iNxn^y Thus 



In particular, 'g\ux{o,t) £ L'^{0,t; L^{U;] 

again, there exists a finite set / C N such 



by finitely many sets UJT' x (0, t 

Let C e Cc°°(int(G);M^''"). Applying Lemma 
that supp(C) ^ Ufce/ ^aT'' ^ (0>^mfc) =: f^- By a partition of unity argument over U, (21) 
holds for g 



4.3 



(b) (ii) 



i): Let t £ (0,T) and U CC be an arbitrary open subset. By Lemma 4.3, we find 



a finite set I C N such that U C IJ^^j and t^fe ^ The claim follows. 



If a relatively open set C is admissible with respect to D we can construct a fine representation 
for H with Lipschitz domains in the following sense. 

Lemma 4.5 (Lipschitz cover of admissible sets) Let H (1 Q be relatively open and admissible 
with respect to D. Then there exists a fine representation {Um} for H such that 

(i) Um is a Lipschitz domain for all m £ N, 
(ii) W~^{dUm n D) > for all m£ N. 

Proof. We will sketch a possible construction for reader's convenience. 

We assume w.l.o.g. that H is path-connected because H can only have at most countably many 
path-connected components and for each component we can apply the construction below. 

Let us choose a reference point xq £ D H with the property 

n''-^ (d{Be{xo) nn)r\D^ > O for ah e>0, (22) 

which is possible since 'H'^^^^Dn H) > 0. The relatively open subset Dm C for m G N is defined 

as 

Dm:=H\Bym{^\H). 

If m is large enough we have xq G Dm since H Q ^1 is relatively open. We define 

D'j^ := {x G Dm I X is path-connected to xq in Dm}- 

Hence, we obtain an e > such that ^^(xo) H 17 C D'^ since D'^ is relatively open in 0. In 
combination with (22), this yields ?^"'~^(9Z)'„ (ID) > 0. Because of CC H, there exists a 
Lipschitz domain Um C O with C Um ^ L[ (e.g. the part of the boundary dUm \ d^l of Um can 
be constructed by polygons such that dUm fulfills the Lipschitz boundary condition). The family 
{Um} satisfies all the desired properties. □ 



Corollary 4.6 Let G C be a shrinking set where G{t) is admissible with respect to D for all 
t G [0,r]. Furthermore, let {tm} C [0,r] be a dense sequence containing T. 

Then, there exists a countable family {UJP'}ke'H of Lipschitz domains UJP' CC G{tm) for each 
m G N such that 

(i) ■W'-^idU]^ r\D)>0 for all m G N, 

(ii) {UJJ^}ken is a fine representation for G{tm) for all m G N, 
(^ii) G = U™=iG(t™) X [0,M. 
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4.1.2 F-limit of the regularized energy 



The construction of the values in (14) satisfying the lower energy bound (16) is based on F- 
convergence techniques which will be introduced below. We refer to |BMR09] for the utilization of 
F-convergence in the context with rate-independent complete damage models. 

Definition 4.7 (F-Iimit of the e-regularized reduced energy) 

Let ■■ X Wl'^in) Moo be fore>0 the (regularized) reduced free energy defined by 



oo else. 



Then, we denote by C; the T -limit of (Bg as e ^ 0+ with respect to the topology in //^(il;M") x 
Wl'^in). Here, W„'^(J^) denotes the space W'^'^{^) with its weak topology. 

Remark 4.8 The existence of the T-limit above is ensured because {G;e} is non-negative and mono- 
tonically decreasing as e — ?■ 0^. Furthermore, ^ is the lower semi-continuous envelope of <to in the 
i?i(Q;M") X Ww'^(rj) topology (see \Bra02^ ). 



To prove properties of the F-limit ^ which are needed in Section 4.3, we will establish explicit 
recovery sequences. The proof relies on a substitution which is introduced in the following. 

Assume that u £ H^{n;W) minimizes J-"e(e(-),z) with Dirichlet data ^ on D. Then, by ex- 
pressing the elastic energy density in terms of its derivative W^^, i.e. = \W% : e, and by 
testing the Euler-Lagrange equation with C, = u — u for a function u E H^{yi;W^) with 5 = ^ on 
D, the elastic energy term in £^ can be rewritten as 



/ W'{e{u),z)dx= [ l{g{z) + e)Ce{u) :e{u)dx. 
Jn Jn 2 



(23) 



Lemma 4.9 For every G H^{Q) and z £ VF^'^(ri) there exists a sequence — )• O"*" such that 
(z — (Je)"'') — )• z) is a recovery sequence for — > 5^ as e — ?• O"*" where 5 is the T-limit of 



given by 

min TM^ + 0,z) ifO<z<l, 
(;eH}y{n;m") 



oo 



else. 



with 



Fe{e,z) := / W'{e,z)dx. 
Jn 



in the H\n;W) x Wi'^{n) topology. 



Proof. The F-limit 5 exists by the same argument as in Definition 4.7 Let {(,£,Zs) — s- {£,,z) be a 
recovery sequence. Since — )• z in C^'°'{Q) due to the compact embedding W^'P{^1) ^ C'^'"($7) 
for some 0<a<l — ^, we can choose a sequence Se — )• 0^ such that (z — Se)^ < z^. Note that 
(z — ^e)"*" G VF^'^(il). Consider the arrangement 



{Z - 5,) + ) - Z,) = dei^, {Z - 5e 
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We observe that ^4^ < because of (note that {z — 6e)~^ < z^) 

Fe{€{i + C), {Z - 5e) + ) < Fe{e{i + C), ^e) 

for all C G Hl,{^]W). Let Ue,Ve G H\){Q.]W) be given by 

Ue= argmin J"e(e(^ + C), ^e), 
CG//i,(n;K") 

Ve = argmin Te{^{ie + C), ^e)- 
CG/^l,(n;R") 



Applying the substitution equation (23) for with u = and for tj^ with u = Ug, we obtain a 
calculation as follows: 

B, = + Ue),Ze) - Fe{e{ie + ^e), ^e) 

+ II ^(^(^e) + e)Ce(u£ + i^s) ||L2(f^) ||e(^ - is) \\L^{n) 

Using ^£ — ^ ^ in H^{Vl), ^ z in W^''P{Q) and the boundednes of -7^e(e(C + Ue),^^) and 7e(e(^e + 
Ue),^^) with respect to e, we end up with lim sup£^o+ -^e ^ 0- Consequently, taking also into 
account that (^e, Ze) — )■ is a recovery sequence, we obtain 

limsup5^e(^, {z - 5e)+) < limsup5'e(^e,2;e) + lim sup + limsup-Be < '^{£,,z). 

e^0+ e-s>0+ e-s>0+ £-5>0+ 



□ 



Corollary 4.10 (i) For every ^ G H^{Q.;W) and z G 



(^iij r/ie recovery sequence {i,{z — 5^)^) — >■ for d constructed in Lemma 4-9 is a 



p 

recovery sequence for — > (S as well. 
(Hi) Let C G H\n;W''), z G W^'P{n) and F C Q be open such that Ipz e W^'P{n). Then 

Proof. 

(i) Let (^e,Ze) —5- ('^, 2;) be a recovery sequence for — > ^. Hence, — ^ in H^{Q;W^) and 
^ z in l^^'^(r2). Applying "liminf£_>Q+ " on each side of the identity 

^s{Ce,Ze)= [ ^\Vz,\Pdx + de{^e,Ze) (24) 

Jn P 
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yields for a subsequence 

(BiC,z) > [ -\Vz\Pdx + ^{tz). 
Jn P 

The "<" - part can be sho wn b y considering a recovery sequence (£, {z — e)^) — s- iS,,z) for 

and applying "liminf£^o+ " ™ (p4|) with {^^, Ze) = (£, (z— e)+) 



r 

— 7- 5^ according to Lemma 



4.9 



on both sides. 

(ii) This follows from (i). 

(iii) Without loss of generality, we assume < z < 1 on $7. Let {z — 5e)^) — t- (£, z) be a recovery 

sequence for — > C; as in (ii). By assumption, \f{z—5£)'^ G W^'P{Q,) and Je)"*" — )• li^z 

in W^^P{a) as e ^ 0+. 

Since fe(e(£ + C), - < £e{e{i + 0,{z - ^e)'^) for all C G Hl{n;W), we obtain 

inf £:,(6(e + C),li^(^-4)+)< inf f,(e(e + C), - 4)+)- 

CeHi,{Q;R") CeHi,{Q;M") 

Therefore, 

Passing to e — )• 0+ yields the claim. □ 



Lemma 4.11 Let £ e M") and z E W'^'P{n) with < z < 1. Furthermore, let u € H^ocii^ > 

0};M") and for every Lipschitz domain U CC {z > 0}, u = on D dU in the sense of traces. 
Then 

^{i,z)<8{e{u),z). 

Proof. Consider an arbitrary e > and define z^ '■= {z — e)^. Since z G C(ri), it holds the compact 
inclusion {z^ > 0} CC {z > 0}. There exists an open set U with Lipschitz boundary such that 
{zgr > 0} C {/ C {z > 0} (e.g. construction of dU \ dO, by polygons such that dU fulfills the 
Lipschitz boundary condition). 

Now, we have u\u £ H^(U;W^) as well as u = £ on dU n D. There exists an extension 
Ue G H^{Q,;W^) with Us\i/ = u\u and = £ on D. The monotonicity of {^e} with respect 
to e implies that (£ is the lower semi-continuous envelope of <B{(,,z) := inf£>o ^£(£, z). in the 
i/i(il;M") X Ww'^(f^)-topology (cf. |Bra02j V By switching the infima, it holds 

^ f infce/fi,(n;R") S (e(C + 0, ^) if < z < 1, 
1 CO else. 

Since u = on {zg > 0}, we get 

£(£,2;)= inf inf liminf (B(£e,Ve) 

(e^S, in m{n;W") ri.^z in VKi'P(n) e-s>0 

< liminf £ (£, z^) < liminf £(€(us),z^) 
e-s>0 e-s>0 



< liminf iS(e(u),Ze) =£(e(u),z). 

e->-0 



□ 
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4.2 Simplified problem 



In the first step of the proof of Theorem 3.10, an existence result of a simphfied problem, where no 
exclusion of material parts are considered, will be shown. The statement we are going to prove in 
this subsection is given as follows. 

Proposition 4.12 (Degenerate limit) Let b G W^''^{0,T;W^'°°{n;W^)) and z° e W^'P{n) with 
<z° <l and {z° > 0} admissible with respect to D be initial-boundary data and let W be given 
by ^ satisfying ^ . Then there exist functions 



z 



G L^iO,T;W''PmnHHO,T;L\n)), u G L^.H^.^^i^oiU > 0});M"), 



eGL2({z>0};M3Tm) 

with e = e{u) in ^£){{z > 0}) such that the properties (ii)-(v) of Definition |3.5| are fulfilled for 
F := {z > 0}. Moreover, (see energy inequality (|14[) ) can be chosen to be it{b^,z^) which 



satisfies (16) by Lemma 4-11 



Remark 4.13 Let us consider the functions e, u and z obtained above in the degenerate limit. We 
do not know that F = {z > 0} equals ^£){{z > 0}) and, if F \ 21d({z > 0}) ^ ill, it is not clear 
whether u can be extended such that e = e{u) also holds in F. On the other hand, we would like 
to stress that (n, z*) with the truncated function z* := zl2tj-,({z>o}) also do not necessarily form a 



weak solution in the sense of Definition 3.5. Because z* viewed as an SBV {Q,T;L {Q))- function 



may have jumps which needs to be accounted for in the energy inequality (|14|). The construction of 
weak solutions will be performed in Section\42 



Let {b^,z^) — )• {b^,z^) with := {z — Se)^ and 6° := 6(0) be a recovery sequence of — > £ 



according to Lemma [4.10 (ii). A modification of the proof of Theorem 4.6 in [HKllj yields the 
following result. 

Theorem 4.14 (e-regularized problem - incomplete damage) Let e > 0. For the given 
initial-boundary data z^ G W^''P{Cl) andb G W^'^ {0,T;W^'°° {^l-jR"^)) there exists a pair = {u£,Zs) 
such that 

(i) Trajectory spaces: 

Ze G L°°(0,r;VFi'P(O))n/7^(0,r;L2(O)), Ue G L°°{{),T;H^{VL;W)). 

(ii) Quasi-static mechanical equilibrium: 

W%{e{u,{t)),Ze{t)):e{Odx = Q (25) 

for a.e. t G (0,T) and for all C, G i?^(r2;M"). Furthermore, Ug = b on the boundary Dt. 
(Hi) Damage variational inequality: 

[ \Vz,{t)r^Vz,{t)-VC + W^MMt)),Ze{t))Cdx> [ {a-(3dtz,{t)-r,{t))Cdx, (26) 
Jn Jn 

Ze{t) > 0, 
dtZeit) < 
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for a.e. t G (0,r) and for all C € W^^P{n) with C < where r^ G L^{Qt) satisfies 

T,{t)i^-z,{t))dx<0 



n 

for a.e. t E (0,T) and for all G W^'^{^1.) with ^ > 0. The initial value is given by 
z^{t = 0) = z° in n. 

(iv) Energy inequality: 

£e{e{Ue{t)),Ze{t))+ [ a\dtZs\ + f3\dtZe\^ d{x, s) 

<£M4),zl)+ f W%{e{u,),z,):e[dth)d{x,s) (27) 
holds for a.e. t G (0,T) where minimizes 

£e{e{-),Ze) in H^{Q.;W) with Dirichlet data 

on D. 

Moreover, r^ in (iv) can be chosen to be 

rs = -XeW,MUe),Ze) (28) 

with Xe £ L'^i^) fulfuHing Xe = ^ on {z^ > 0} and < Xe ^ 1 on {z^ = 0}. 

We consider a sequence {smImsn ^ (0, 1) with em — )• 0"'' as M — )• oo and for every M G N 



a weak solution (ug^,^, z^j^^) of the incomplete damage problem according to Theorem 4.14 The 
index M is omitted in the following. We agree that := e{us) denotes the strain of the regularized 
system. Our further analysis makes also use of the truncated strain eg (the strain on the not 
completely damaged parts of Cl) given by 

Ce := eet{z^yo}. 

We proceed by deriving suitable a-priori estimates for the incomplete damage problem with respect 
to e. 

Lemma 4.15 (A-priori estimates) There exists a C > independent of e such that 

(i) l|e.||L2(n^.Knxn) < C, (iii) WdtZeWmnr) ^ 

(ii) SUPte[o,T] \\Ze{t)\\w-^,P(n) < C, H W^i^^e, 2e)||L-{0,T;Li(n)) < C. 



Proof. Applying Gronwall's lemma to the energy estimate (27) and noticing the boundedness of 
<5e(e(u^), z^) with respect to e G (0, 1) show (iii) and 

£,{ee{t),Ze{t))<C (29) 

for a.e. t G (0,T) and all e G (0,1) (cf. |HK11) ) and in particular (iv). Taking the restriction 
< < 1 into account, property ([29]) gives rise to ||-Ze||Loo(o,r;Wi'P(r2)) ^ C'. Together with the 
control of the time-derivative (iii), we obtain boundedness of ||-Ze(t)||vFi'P(fi) ^ ^ every t G [0,T] 
and e G (0, 1). Hence, (ii) is proven. 
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It remains to show (i). To proceed, we test inequality (26) with = — 1 and integrate from 
t = to t = T: 



/ W^,{ee,Ze) +red{x,t) < a- Pdtzed{x,t). 
Applying ([6]), ([28]) and ([30]), yield 

J{z,>0} 

- I 7;9'{z£)Ce£ : e£d{x,t) 

J{ze>0} ^ 

= [ W%{ee,ze)d{x,t)- I W';,{ee,ze)d{x,t) 



< / W^,{ee, Ze) d(x, t) - / XeW^,{ee, d(x, t) 



WUee,Ze)+reit)d{x,t) 



< / a - PdtZsd{x,t). 



This and the boundedness oi a — /3 dtZird{x,t) with respect to e shows (i). 



(30) 



□ 



Lemma 4.16 (Converging subsequences) There exists functions 

e e L2(!^r;M"^"), z E L°°(0, T; W^'P{n)) n H\0, T; L'^{n)), 

where z is monotonically decreasing with respect to t, i.e. dtz < 0, and a subsequence (we omit the 
index) such that for e — 0+ 



(i) z.^zin F1(0,T;L2(S7)), 
z.^zin LP{0,T;W^^P{n)), 
Ze{t) z{t) in W^^P{Q.), 



(ii) Ce^ein ^^(Ot; M""""), 

W%{ee,ze) W,e{e,z) in L'^{{z > 0};M"^"), 
Wl{ee,Ze) ^ in L'^{{z = 0};M"^"). 



— )• z in Q 



T, 



Proof. The a- priori estimates from Lemma 4.15 and classical compactness theorems as well as 
compactness theorems from J.-L. Lions and T. Aubin yield [ Sim86j 



Ze^zm L°^(0,r;TFi'P(J7)), 
Ze^ z\n H\0,T;L'^{n)), 
Ze ^ z in LP{Qt) 



W^,iee,Ze)^We in L2(J7r; M"X"), 



as e — )■ 0+ for a subsequence and appropriate functions We, e and z. 

Proving the strong convergence of Vz^ in LP(Ot;]R") does not substantially differ from the 
proof presented in [HKllj . It is essentially based on the elementary inequality 

Cnc\x - y\P < {{\xr^X - \y\P-^y),x- y) , 
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where (•, •) denotes the standard euchdean scalar product and on an approxhnation scheme {Q} ^ 
LP{0,T;W^^P{Q)) with Q > and 



Ce ^ z in LP{0, T; W^'P{n)) as e ^ 0+, 
0<Ce<Zs a.e. in Qt for all e G (0, 1). 



(31a) 
(31b) 



Using the above properties, we obtain the estimate: 

Cue/" \VZe-Vz\Pd{x,t)< [ {\Vz,J'^VZe-\Vz\P-'^Vz) -ViZe- Z)d{x,t) 

= I |Vz,|P-2Vze-V(ze-C.)d(x,t) 
Jq.t 

^ ' 

+ f \Vze\P''^Vze ■ V{Ce - z) - \V z\P-^V z ■ V {ze - z) d(a;, t) . 



The weak convergence property of {Vz^} in LP[Q.t) and (31a) show — )• as e — )• 0+. Property 
(26) tested with C,{t) = Ce{t) — Ze{t) and integration from t = to t = T yields 



WUe{Ue),Ze){Ce-Ze)dix,t) + 



-a + /3{dtZsit))){Cs- Ze)d{x,t) 



<0 by ^ and |3Tb]l 



>0 as e->0+ by ||31a| 



Here, we have used r^^ = on Qt (see (|28[)). Therefore, (i) is also shown. 
To prove (ii), we define to be {z^ > 0} n {z > 0}. Consequently, we get 



WUee,Ze)lN,=W'iee,Z,)l 



and the convergence 



1{2>0} in ^1 



(32) 



(33) 



for e — 7- O"*" by using — t- z in Q,t- Calculating the weak L^{Qt]^ 



-limits in ( 32 ) for e — )■ 0"'" 



on both sides by using the already proven convergence properties, we obtain W^e(e,z) = We- The 
remaining convergence property in (ii) follow from Lemma 4.15 (iv). □ 

We now introduce the shrinking set F C by defining 

F{t) := {zit) > 0} 



for all t G [0,T]. This is a well-defined object since F C fij^ is relatively open by Theorem 
well as F(s) C F(t) for all < t < s < T by the monotone decrease of z{x, ■). 



A.2 



as 



Corollary 4.17 Let t G [0,r] and U CC F{t) be an open subset. Then U C {ze{s) > 0} for all 
s G [0, t] provided that e > is sufficiently small. More precisely, there exist < £o,r] < 1 such that 

Zeis) > rj in U 



for all s G [0, t] and for all < e < eq. 
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Proof. By assumption, we obtain the property dist{U, {z{t) = 0}) > 0. Therefore, and by 
z{t) E C{Q), we find an r/ > such t hat z{t) > 2rj in U . By exploiting the convergence Ze{t) — )• z{t) 
in as e — )■ 0"*" by Lemma 4.16 (b) and the compact embedding ^ C(r2), there exists 



an eo > such that Ze{t) > rj on U for all < e < eo- Finally, the claim follows from the fact that 
Ze is monotonically decreasing with respect to t. □ 



Lemma 4.18 There exists a function u G L^ff^ j^^(2lz)(F); M") such that 
(i) e{u) = e a.e. in^£,{F), 
(a) u = b on the boundary Dx H ^£){F). 



Proof. Let {f/™} and {tm} be sequences satisfying the properties of Corollary 4.6 applied to %d{F)- 
We get for each fixed A;, m € N 



uj^x [o,M c{z, >o} 



for all < e ^ 1 due to Corollary 4.17 Inclusion (34) implies 

e{Ue) = Ce 



(34) 



(35) 



a.e. in f7™ x (0,tm)- Korn's inequality applied on the Lipschitz domain C/™ yields (note that 

n''-\du]^ n D) > 0) 



l'"e|lL2(0,t,„;Hi(C/™;R")) 



< 2 



\Ue{t) - ?>(t)||ii([/™.iRn) + ||fc(t)||^i(C7- 



dt 



<C 1 + 



k (*) ) 1 1 L2 ((7™ ;IR" X " ) 



dt 



<C 1 + 



leeim 



dt 



with a constant C = C{Ujr', b) > 0. Together with the boundedness of Ce in L'^{Q,t] I 
find a subsequence e — )• 0"*" and a function ■u^'^'™') e L?'{Q,tm] H^{Uf^;W^)) such that 



we can 



(36) 



Thus e(u'^'^'™')) = e in x (0,tm) because of (35) and the weak convergence property of e^. 
For each k,m £ N, we can apply the argumention above. Therefore, by successively choosing 
subsequences and by applying a diagonalization argument, we obtain a subsequence e — )• 0"'" such 



that (|36|) holds for all A;, m e N. 

n(*=2'"^2) a.e. on C/^^ x (0,t^J n [/^^ x (0,t^J for all fci, fca, mi, ms G N, 



Since u^^^'"^^^ 

we obtain an n : ^d{F) M" such that n|c/-x(o,t™) ^^(0, t^; i^H^^r; J^")) ^h m G N. 
Proposition 4.4 (a) yields u £ L^H^^^^{F;W^) and the symmetric gradient e(n) coincides with e. 
Therefore, (i) is shown. 

Furthermore, for every A;, m G N, we have u(t) = 6(t) on c?t/^ n D in the sense of traces for a.e. 
t G [0,tm]- By Proposition |44](b), (ii) follows. □ 



We are now able to prove Proposition 4.12 
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Proof of Proposition \4-lS\ Lemma 4.16 and Lemma 4.18 give the desired regularity proper 



4.12 



ties of the functions (e, u, z) in Proposition 
property e = e{u) in %d{F) follows from Lemma 4.18 



Here, we set e := €\f e L^(F;M"^"). The 



In the following, we are going to prove that properties (ii)-(v) of Definition 3.5 are satisfied. 



(ii) Lemma 



4.16 



ii) allows us to pass to e — )• O"*" in (25) integrated from t = to t = T. 



Therefore, equation (11) holds for a.e. t G (0,T) and all Q G M"). Moreover, the 

boundary condition u = 6 on Dt H %o{F) is satisfied. Definition 3.1 immediately implies 



(iii) We first show (fT2|. Let C e -L°°(0, T; VF^'P(O)) with C < 0. The variational inequality ^ 



and the representation for (28) imply 



0< / \V Zer^Vze-VC + {-a + /3dtZe)C dix, t) + [ WUee,Ze)Cdix,t). (37) 



In addition. 



/ W^,{ee,z,)Cd{x,t)< I W%{ee,z,)Cd{x,t) 

J{ze>0} J Fn{ze>0} 



4.16 



Lemma 



proof of Lemma 4.16) yield 



a lower semi-continuity argument and l{z^>o}n{z=o} ~^ l{z=o} ^-S- in f^T (see 



limsup /" W%{e,,Ze)Cd{x,t) < [ Wl{e, z)C d{x,t). 

e->0+ J{z^>0} ' JF 



Therefore, applying "limsup£_^Q+ " on both sides of (37), using the above estimate and Lemma 



4.16 yield 



/ \Vz\'P-'^Vz-VC + W^^{e,z)Cdx> [ {a - l3dtz)C dx. 
JF ' Jn 



(38) 



The properties dtz < and z > a.e. in J7t follow from Lemma 4.16 by taking Stz^ < and 
Zs > a.e. in into account. 



(iv) The jump condition (13) in (iv) of Definition 3.5 holds trivially since we have the regularity 
z E L°°{0,T;W^'P{n))nH^{0,T;L'^{n)). 



(v) To complete the proof, we need to show the energy estimates ( 14 ). Since {b^, z^} is a recovery 
sequence, we get £s{^{u^), z^) — >• (B{b^,z^) as e — 0"*". Now, applying " lim sup£_^o+ " on both 



sides in (27) and using the convergence properties in Lemma 4.16 as well as lower semi 



continuity arguments yield 



Jo Jf{s) 



e,z) : e{dtb)d{x,s) 



> limsup i£eiee{t),Zeit))+ / a\dtZe\+ P\dtZe\'^ d{x,s)] 

E^0+ \ JUt / 
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>limsup / W'{esit),Zeit))dx+ [ -\Vz{t)\P dx 
£-s>o+ Jn Jn P 



+ / a\dtz\+ f3\dtz\'dix,s). 
Jnt 



Indeed, for an arbitrary t G (0, T), we derive by Fatou's lemma and Lemma 4.16 
/ I limsup / W^{es{s),Ze{s))dx] ds > limsup / {ce, z^) d{x, s) 

> liminf / {g{zs) + e)Ce£ : Ce d(x, s) 



> / W{e,z)d{x,s). 
Jf 



We have used the weak convergence property 

y/gize) +eee \/ 9{z) e in L^(0t;1 



as e — )• 0"*". To the end, (40) imphes 



hmsup / W^{ee{t),ze{t))dx> / W {e{t) , z{t)) dx 
£-s>o+ Jn JF(t) 



(39) 



(40) 



for a.e. t G (0,T). Combining it with (39), estimate (14) is shown. 



□ 



4.3 Existence of weak solutions 

By using the achievements in the previous section and Zorn's lemma, we will prove the main results. 
Theorem 3.10| and Theorem |3. 11 To proceed, let > be fixed and V be the set 

T' := { (T, e, u, z, F) I < T < T and (e, n, z, F) is an approximate weak solution on 
[0, T] with fineness r} according to Definition 3.8 }. 

We introduce a partial ordering < on P by 



(Ti, ei. Ml, zi, Fi) < (r2, 62,^2,^:2,-^2) <^ Ti<T2,e2 

^2l[0,fi] = 

The next two lemma prove the assumptions for Zorn's lemma. 
Lemma 4.19 P / 0. 



[0,Ti] 



ei, U2 



[0,Ti] 



Ui, 



Proof. Let (e, u, z) be the tuple from Proposition 4.12 to the initial-boundary data (z , b). If there 



exists an e > such that Jz* H [0,e] = with z*{t) := z{t)l(^^(^^z-{t)>o}) then {e, e,u, z, F) G V. 
Otherwise, we find G Cz*- We claim 



({z° > 0} \ 2tD({z(t) > 0})) ^ as t ^ 0^ 



(41) 
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We consider the non-trivial case z ^ 0. Let X G {z° > 0} n Since {z^ > 0} C fij^ is relatively 
open and admissible with respect to D, there exists a Lipschitz domain U CC {z^ > 0} with x E U 
such that 'H"~^((9C/ HD) > by Lemma 4.5 Because of Theorem A. 2 , z G C{Qt) and, consequently, 



there exists a t > such that U CC {-^(s) > 0} for all < s < In particular, x £ 21d({-z(s) > 0}) 
for all < s < t. This proves (41). Finally, choose e > so small such that e < tj and (note the 
monotonicity of z with respect to t) 

{{z{t) > 0} \ ^D{{z{t) > 0})) < {{z' > 0} \ ^D{{z{t) > 0})) < V 

for all < t < e. We have proved that {e,u,z) on F := {z > 0} is an approximate weak solution 
with fineness ry on the time interval [0, e], i.e. {e, e, u, z, F) £ V. □ 



Lemma 4.20 Every totally ordered subset ofV has an upper bound. 

Proof. Let 7^ C P be a totally ordered subset. We denote with [0, Tr] the corresponding time 
interval of an element R £ TZ. Let us select a sequence {Tq, eg,ug, zq, -Pelggj-g ^ ^) with Tg-^ < Tg^ 

for 6*2 < Oi and linig^o+ ^g = supg^-;^ Tq—.T. 

Let t £ (0, T). There exists a 9 £ (0, 1) with Tg > t and we define 



(e(t),n(t),z(t),F(t)) := {eg{t),ugit), zg{t), Fg{t)). 

By construction, the functions {e,u,z) satisfy the properties (ii)-(v) of Definition 
remains to show that {e{t),u{t), z{t)) are in the trajectory spaces as in Definition 
F satisfies Definition |3.8| (ii). 

The energy estimate for {eg,ug, zg) implies 



3.5 



3.8 



on [o,r]. It 

(i) and that 



£{e{t),z{t)) + 
<4 + 



a\dfz\+ ^\dfz\^d{x,s) 

JF{s) 

W,e{e,z) : eidtb)dix,s) 



JF(s) 



(42) 



for a.e. t £ (0,T). Gronwall's lemma yields boundedness of the left hand side of (42) with respect 
to a.e. t £ (0,f). 

We immediately get 



z £ L°°(0, f ; W^'Pin)) n SBV\0, f; L^{VL)), 



(43) 



Variational inequality (12) tested with (, = —1 shows 

W^-^{e{t),z{t))dx < adx 
F{t) ' Jn 



f3dfz{t) dx 



for a.e. t £ {0,T). This implies 



e £ L\F- : 



(44) 



We know that u\ux{o,t) e L'^{0,t; H^{U;W)) for ah t £ (0,f ) and all open subsets U CC 2tr»(F(t)). 
Let {Uk} be a Lipschitz cover of the admissible set 



F{T) ■.= ^d{{z'{T)>0}) 
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according to Lemma |4.5| (in particular, Definition 3.8 (ii) is fulfilled). For each /c G N, we apply 
Korn's inequality and get for all t £ (0, T) 

\W - ^\L^{0,t■,H^U^:■,R")) ^ C'lk(^)llL2(0,t;L2(c/j^;Rn)), 

where C > depends on the domain but not on the time t. Thus ^I^^^^q £ L'^{0, T; H^{Uk; M")). 
In conclusion, 



(45) 



Therefore, property (i) of Definition |3.8| follows by (43)- (45). We end up with 
{T, e, li, z, F} G V satisfying {Tq, eg, ug, zg, Fg} < {T, e, n, z, F} for all 9 G (0, 1). □ 

Weak solutions exhibit the following concatenation property. 
Lemma 4.21 Let ti < t2 < be real numbers. Suppose that 

q := {e,u,z,F) is an approximate weak solution on [ti,t2], 
q := {e,u,z,F) is an approximate weak solution on [t2,t3] 

with'^l^ = G;(6(t2), ^^(^2)) (ihe value for q in Definition 3.5). 

Furthermore, suppose the compatibility condition z'^{t2) = z^ {t2)'^<^^(^{z~ (t2)>0}) ^'^^ Dirichlet 
boundary data b G W^'^{ti,t3;W^'°°{Q;M.^)). Then, we obtain that q := {e,u,z,F) defined as 
l\[ti,t2) ■~ ^ ^^'^ '?l[t2,i3] ? '^'^ approximate weak solution on [tijts]. 



Proof. Applying "lim^_j_^- ess inf^g^g j^)" on both sides of the energy estimate (14) for {e,u,z,F) 
yields 

lim essinf ^(e(r),z(r)) + / ' /" a|5fz| + /3|9f zj^ d(x, s) + liminf V Jr 



t2 



'h JF{s) 

This estimate can be rewritten as 

rt2 



W4e,z) : €{dtb)d{x,s) 



(t{bit2),Z+{t2)) + 



ti JF(s) 



a\dfz\+ ^\dfz\^d{x,s) 



+ liminf V + lim ess inf £:(e(r), z(t)) - e(6(t2), 2+(t2)) 



TeJzn{t2, 



t2 



W,eie,z) : eidtb)d{x,s). 



(46) 



ti JF{s) 



In the following, we show that we may choose the value (B{b{t2) , z~^ {t2)) for e^^. By the property 
(i) of Definition [3^ we get z'is) ^"(ta) in W'^^P{n) and b{s) b{t2) in H^i'°^(17; M") as s ^ t^ . 
In particular, by using Lemma |C.l and the monotone decrease of z~ with respect to t, 

Z~{s)l<^n{{z-{s)>0}) ^"(*2)ln,g,^^_^^)2la({^-(r)>0}) =: X 
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in W^'^{^) as s — ?• • the definition of Xi tlie inclusion 



2tD({^"(t2) >0}) C Pi ^d{{z-{t)>0}) 

TG(tl,i2) 

and tlie compatibility condition, we find z^(t2) = xlao({z-(i2)>o})- 

Thus, applying Lemma 4.11, lower semi-continuity of the F-limit ^ and Corollary 4.10 (iii), we 
obtain 

lim ess inf £^(e(r), z(r)) = lim essinf <f (e(T), z^(r)) 



s-5>t2 -re(s,i2) 



This leads to 



> lim essinf £:(e(n(T)),z (r)l<ao({2-(r)>o})' 

S-S>t~ T6(s,t2) 

> lim essinf £(6(r),z-(r)l2t^(|^_(^)>o})) 

s^t2 •re(s,t2) 

>mt2),x) 



0< V < lim essinf ,f(e(r),z(T)) - e(6(t2),2+(t2)) 

seJzn(ti,i2] 2 



-l-liminf ^ J'^, 



2 TeJ2n(ti, 



where the second '<' becomes an '=' if f2 G Jz- Consequently, (46) becomes 



(t{b{t2),Z+it2)) + 



t2 



ti JF{s) 



a\dfz\+ I3\dfz\^d{x,s) 



t2 



seJzn{ti,t2] 
W,e{e,z) : eidtb)d{x,s). 



(47) 



The energy inequality (14) for {e,u,z,F) (taking "e^ = <S.{b{t2),z'^{t2)) into account) can be ex- 
pressed as 

£{e{t),z{t))+ [' [ a\d!z\+P\dfz\^d{x,s)+ V Js 

Jtn JF(s) _ rrf. .1 



'i2 JFis) 

< e(6(t2),Z+(i2)) 



t2 -'i^(s) 



s&j^r\{t2,t] 
W4e,z) : e(9t6)d(x,s) 



(48) 



for a.e. t E (t25i3)- Adding (47) and (48) shows that the energy estimate for {e,u,z,F) also holds 
for a.e. t G (i2,t3)- It is now easy to verify that {e,u,z,F) is a approximate weak solution on the 
time interval [^1,^3] according to Definition 3.5 □ 



Proof of Theorem \3.10 . By Zorn's lemma, we deduce the existence of a maximal element 
R = {T,e,u,z,F) in V. In particular, a maximal element satisfies the properties in Theorem 



3.10 



on the interval [0,T]. We deduce T = T. Otherwise, we get another approximate weak solu- 
tion (e, u, z, F) on [T, T+e] for an e > with initial datum z~ {T)lf^^^~_ ^j;-^^^-^ (which is an element 
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of W^'P{n) by Lemma 
Lemma 



4.21 



with e± = €{h{T),z{T)) if r G J,. By 



C.l ) as in the proof of Lemma 4.19 

(e, u, z, F) and (e, u, z, F) can be concatenate to an approximate weak solution on 
[0, T + e] which is a contradiction. □ 



Proof of Theorem 3.11\ Here, let us consider the set V given by 

V := {(r,u,z) |0 < T < T and (n, z) is a weak solution on 
[0, T] according to Definition 

with an ordering < as above (except the conditions 62]^^;^] = 



3.5 



ei and i^2|[o,fi] 



which are 



not needed here). Proposition 4.12 shows 7^ 7^ by noticing z G C{Qt) (see Theorem A. 2) and 
< T] < z^. The property that every totally ordered subset of V has an upper bound can be shown 
as in Lemma 14.201 A maximal element satisfies the claim. □ 



A Embedding Theorem 

The embedding theorem |A.2| in this appendix is a special version of a more general compactness 
result in |Sim86| Corollary 5]. However, we would like to present a different (short) proof which 
requires the following generalized version of Poincare's inequality. 

Theorem A.l (Generalized Poincare inequality |Alt99^ Section 6.15]) Let Q C be a 

bounded Lipschitz domain and M C W^'^{Q;W^) non-empty, convex and closed with 1 < p < 00. 
Furthermore, M satisfies the property 

u e M, a>0 =^ au e M. 

Then the following statements are equivalent: 

(i) There exists a uq £ M and a constant Cq > such that for all ^ G M*" 

UQ + ^GM =^ \^\<Co. 

(a) There exists a constant C > such that for all u £ M 

LP {n;R"^) ^ C*!! V^IIlp(Q;K"X")- 

Theorem A. 2 Let Q M" be a bounded Lipschitz domain and p > n. Then 

L°°(0, T; W^^P{n)) n H\0, T; L^{VL)) C C{Tl^). 

Proof. Let z G L°°(0, T; W^'P(Q)) n iJi(0, T; ^^(il)). We can choose a representant such that z G 
C([0, T]; L2(0)) and z{t) G W^^P{n) for all i^[0, T]. By employing the embedding W^'P{^1) C 0(^1) 
(note that p > n), we obtain a function z : f^T — )• M such that 

z G C{[0,T];L^{n)) and z{t) G C{U) for ah t G [0,r]. (49) 

Let ( 

Xmit-rri) S be arbitrary with {xuntm) — ^ {x,t) in Q,x as ttt. — > 00. We have 

|2;(x, t) z(^X-iYi} im)| ^ |^(''^) ^) ^(.■^1 tm)\ ~l~ [-^(s^, ^m) ^(^rrt) ^m)| ■ 
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Assume that 7^ as m — )• cxd. Then there exists a subsequence of {^m} (also denoted by {Am}) 
such that hmm->oo Am > 0. Using this subsequence, it holds z{-,tm) — ^ z{-,t) in L'^{Q) due to (49). 
We obtain again a subsequence (we omit the additional subscript) such that z{y,tm) z{y,t) as 
m — )■ 00 for a.e. y G il. Therefore, we can choose ym ^ x in Q, such that \z{ym, t) — z{ym, tm)\ — ^ 
as m — )• 00. It follows 

\z{x,t) - z{x,tm)\ < \zix,t) - z{ym,t)\ + \z{ym,t) - z{ym, tm)\ + \ziym, tm) " z{x,tm)\ 
^ V ' ^ V ' ^ V ' 

41 A2 AS 

-'^m ■'^m 



The continuity of z{-,t) due to (49) implies A^ — )• as m — )• 00. A,^ converges to by the 
construction of {ym}- To treat the term A^, we apply the Poincare inequality in Theorem A.l to 
M := {u G W'^^P{Bi{qo)) \ u{qo) = 0} and obtain 

< C\\Vg\\ (50) 

for all g G VF^'^(i?i(go)) with g(go) = 0, where go G I^" and C > is independent of g and qo. Note 



that, due to 5 G W^''P{Bi{qQ)) C C(i?i(go)), 9 is pointwise defined. By utilizing ( |50| ) and using a 
scaling argument, we gain a C > such that for all e > and all g G l^^'P(i?£(go)) with g{qQ) = 
follows 

p — n 

MciBAm)) = II^^(^-)IIc(bTM) - ^ll5(^-)llw'i.p(iJi(<?o)) < C'||eV5(e-)llLf(Bi(go)) = C^^" l|V5llLP(B,(go))- 

By setting gm{-) '■= z{ym,tm) — z{-,tm) and Sm : = 2 1 — 2; | , we can estimate A^ in the following 
way (note that gmiUm) = 0): 



p — n 



Since z G L°°(0, T; T^1'P(J7)) n if^O^ ^^(^^))i W"^ 9m\\Lv{B,^(y^)) is bounded with respect to m. 
In conclusion, A^ — )• as m — )• 00. Hence, we end up with a contradiction. Therefore, Am — )• as 
m — )• 00. 

The convergence — >• as m — )■ 00 can be shown as for A^ — >• 0. □ 



B Chain-rule for vector-valued functions of bounded variation 

Theorem B.l (BV-chain rule |MV87| ) Let I (1 M. be an interval, X he a real reflexive Ba- 
nach space, f G BViodl', X) with df = f'^fx for a non-negative Radon measure fi on I and 
f'^ G L^^^I, X). Moreover, let E : X ^ be continuously Frechet-differentiable. Then 
Eo f £ BViocil',^) and d{E o /) admits as density relative to ^ the function t 1— )• {9{t), f'^{t)) , 
where 6:1^ X* is defined as 

e{t) := /' di?((l - r)f{t-) + r/(t+)) dr. 
Jo 

Corollary B.2 Suppose f G SBV{0,T; X) and E : X ^M. is continuously Frechet-differentiable. 
Then Eo f e SBV{0, T) and for all < a < b < T: 

d{E o /)((a, b]) = f\dE{f{s)),f'{s)) ds+ (^(/(«^)) - E{f{s-))) . 

seJ/n(a,b] 
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Proof. We apply Theorem B.l By assumption, we obtain the decomposition d/ = /^/x with 
/i = C^+n^lJf and f'^{t) = T{t) + f{t+) - f{t-) for all t e (0,r). Applying Theorem yields 



d{Eof){{a,b])= / (0(s),/;(s))dM5) 

J{a,b] 



= [ {eis)J'{s))dC\s)+ {0{s)J{s~')-f{s-)) 
•^("•''1 teJfn{a,b] 

Since /(s^) = /(s~) = /(s) for a.e. s G (a, 6], the first term on the right hand side becomes 



{eis),ns)) dC\s) = ( / dE{{l - r)f{s~) + r/(s+)) dr, /'(s)) dC\s) 

(a,b] J{a,b] Wo ' 

= [ {dE{fis),ns))ds, 

J{a,b] 

where ds := dC^{s). Furthermore, by the classical chain rule, 

{eis)Jis+)-f{s-))= Yl {£dEii^-r)f{s-)+rfis+))dr,f{s+)-fis-)) 



seJ/n(a,6] seJfn{a,b] 

= E 

se J/n(a,fe] 



Y f {^m^ - r)f{s-) + rf{s+)), f{s+) - f{s-)) dr 

Jfn(a,b] ''^ 

Y r^i?((l-r)/(s-) + r/(.+ ))dr 



se J/n(a,fe] 

Y {E{f{s^))-E{fis-))). 

s£jfn(a,b] 



□ 



C Truncation property for Sobolev functions 

Lemma C.l Let D,Q C M" be open sets and p > n. Furthermore, assume that a function 
f G T^^'P(O) fulfills / = on dD \ dQ (f is here considered as a continuous function due to the 
embedding W^^P{n) ^ C{U)). Then fin G W^^P{n). 

Proof. We can reduce the problem to one space dimension by using the following slicing result from 
[AFPOOl Proposition 3.105] for functions u £ LP(Q): 

u G W^^P{n) ^ Vi/ G S"-^ : < G W^^PiCf^) for C^'^-a.e. x e 

and / / |V<|Pdtdy < oo, (51) 

where ^li, is the orthogonal projection of 0, to the hyperplane orthogonal to and := {t G 
M I X + tz^ G ri} as well as ti^(t) := u{x + ti^). 
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Applying this result to /, we obtain G W^'P{Q'^) for £""^-a.e. x e Q„ and all u G S"~^ 
Moreover, slices for the function g := fl£) are given by the equation 



9x 



The function is absolutely continuous. We claim that this is also the case for g'^. To proceed, 
let e > be an arbitrary real. Then, we get some constant 5 > such that 



(afc, k G /, with < are finitely many disjoint intervals of 17^ with \ cLk — bk\ < S 

k&I 

=^ J2\rAak) - rAbk)\ < e. (52) 



The property (52) is also satisfied for g^. Indeed, let {ak,bk), k G I, with < bk be finitely many 
disjoint intervals of Q'^ with X^fce/ I'^fc ~ ^^1 ^ ^- define the values ak and 6a; in the following 
way: 



{ak,bk) 



{ak,bk) a ak,bk € or ak,bk ^ D^, 

{z, bk) for an arbitrary fixed z £ [ak, bk] H dD^ if ak and G -D^, 
^ (afc, z) for an arbitrary fixed z £ [ak, bk] H dD^ if ak G -DJJ and Z?J!. 



We conclude Xlfce/ l«fc - bfcl < Efce/ l^fc -bk[ <S and therefore Xlfce/ l/x («fc) - /x (&fc)| < e by ( |52| ). 
Taking 

fce/ fce/ fee/ 

into account, shows that g^ is absolutely continuous and we find g^ G W^'P{^';^). 

Moreover, /^^ iV^^l^dtdy = /^^ |V/^^|Pdtdy < oo. Applying ^ yields g G VFI'P(O). 

□ 
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